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The lightest neutraUno, as the dark matter candidate, can be gravitationally cap- 
tured by the Sun. In this paper, we studied the high energy neutrino signals from 
solar neutralino annihilations in the core of the Sun in the anomaly mediated su- 
persymmetry (SUSY) breaking (AMSB) model. Based on the event-by-event monte 
carlo simulation code WimpSim, we studied the detailed energy and angular spec- 
trum of the final muons at large neutrino telescope IceCube. More precisely we 
simulated the processes since the production of neutrino via neutralino annihilation 
in the core of the Sun, neutrino propagation from the Sun to the Earth, as well as 
the converting processes from neutrino to muon. Our results showed that in the 
AMSB model it is possible to observe the energetic muons at IceCube, provided that 
the lightest neutralio has relatively large higgsino component, as a rule of thumb 
+ > 4% or equivalently cfsd > 10~^p6. Especially, for our favorable param- 
eters the signal annual events can reach 102 and the statistical significance can reach 
more than 20. We pointed out that the energy spectrum of muons may be used to 
distinguish among the AMSB model and other SUSY breaking scenarios. 

PACS numbers: 95.35.+d, 12.60.Jv, 13.15.+g, 95.55.Vj 



I. INTRODUCTION 

Over the past several years, the picture that our universe is not baryon dominant is well 
established. Based on the observations of Wilkinson Microwave Anisotropy Probe (WAMP) 
and other experiments, the relic dark matter density is fixed to be 0.085 < Q^Mh"^ < 0.119 



within 2a uncertainty [l|. In order to account for the dark matter, physics beyond the 
standard model (SM) is usually required. Among these new physics the weakly interacting 



2 



massive particle (WIMP) is one of the popular candidate for the dark matter. One possible 



WIMP is the 
conservation 



ightest supersymmetric (SUSY) particle (LSP) in SUSY theory with R-parity 
^. Usually there are three possible LSPs, namely sneutrino, gravitino and 



neutralino. The lightest sneutrino has been largely ruled out by direct searches 



3| , and the 



Cavitino is difficult to be detected due to its very weak interactions with ordinary matter 
]. Thus the lightest neutralino, which is combination of gaugino and higgsino, is the most 



well-motivated candidate for dark matter 

SUSY breaking mechanism plays an important role to investigate the nature of dark 
matter. Firstly, in order to be consistent with experimental observations, SUSY must be 
broken at weak scale. As a consequence, the mass of dark matter is partially determined 
by the breaking terms. Secondly, we need higher scale SUSY breaking mechanism to reduce 
the number of free parameters at low energy. As an example there are more than 100 
parameters in minimal supersymmetric standard model(MSSM). Thirdly, the properties of 
the LSP depend on the details of SUSY breaking mechanism. For example in minimal 
supergravity (mSUGRA) model, the main component of the lightest neutralino can be the 
bino in 'bulk region' while in 'focus point region', the lightest neutralino can have larger 
higgsino content sl, Q, [s, [ol, [l^. Here 'bulk region' and 'focus point region' are determined 
by different SUSY breaking parameters. Thus it is interesting to explore (i) the predictions 
for different SUSY breaking mechanisms; and (ii) whether the detection of dark matter can 
distinguish among the different SUSY breaking mechanisms. 

While in literature the mSUGRA model was widely investigated, in this paper we are 
interested in the properties of neutralino (as the dark matter par ticle) and its detection in 
the scenario of the anomaly mediated SUSY breaking (AMSB)|lll.[l^. Many predictions in 
AMSB models are different from those in usual mSUGRA models and other SUSY breaking 
scenarios. One feature of AMSB model is the special gaugino mass relation calculated by 
(3 function. Mi : M2 : IM3I = 2.8 : 1 : 7.1 which implies that the lightest chargino 
pair and neutralino comprise a nearly mass-degenerate triplet of winos over most of the 
MSSM parameter space. Moreover the lightest neutralino is mostly wino-like. Such kind 
of neutralino usually has larger cross section compared to that in mSUGRA model. As a 

nn 

consequence, if it is the dark matter, it should be produced in non-thermal processes |13I.I14I|. 
or should be very heavy, generally larger than 2 TeV 

The methods to detect dark matter can be classified into direct and indirect ones. The 
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former methods detect dark matter by measuring the recoil of nucleus. The latter ones 
detect the final stable particles, including neutrinos, antiprotons, positrons, antinuclei and 
photons, which are produced by dark matter annihilation. In this paper we will focus on the 



observation by neutrino telescopes which detect t 
annihilations (for some recent works, see 
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le high ener gy n eutrinos from dark matter 
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25|, 126|, 1271). The 



non-relativistic dark matter would be trapped by Sun or Earth via gravitational force. The 
captured dark matter can annihilate into SM particles, including neutrinos, in the center 
of Sun or Earth 
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30|. The resulting neutrinos interact with usual matter feebly, thus 



it is possible for us to detect them after a long propagation. Many neutrino telescopes for 



high energy neutrino have been built or are under construction such as IceCube 
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Super-Kamiokande 33ll, AMANDA 
on the neutrino detection at IceCube. 
pole and will be completed in 2010 [31 



32|, 



34| and ANTARES [35|] etc. In this paper we will focus 
he IceCube detector now is being built in the south 



32| . The scale of the IceCube is in square kilometers 



3l| 



and it can detect neutrinos with threshold energy of 50 GeV 

In this paper, we will discuss neutrino signals from neutralino annihilations in the Sun in 
AMSB scenario. The final muon events rate at the IceCube detector can be determined by 



36|: 



events 



where Tann is neutralino's annihilation rate in the Sun, A^ff is detector's actual effective 
detected area, Bvi is the branching fraction of i-th annihilation channel, and $^ is the fi flux 
over detector threshold per unit area and parent pair. Our study shows that for the Sun 
there exit the parameter space which can induce the large neutralino capture rate and thus 
the large neutralino annihilation rate Tann- In this case, most neutrinos come from W decay 
and the Ws are pair produced via neutralino annihilations. Such neutralino is wino-like. 



Recently the authors of Ref . 



26l ] pointed out that if the spin effects of neutralino annihilation 



are included, the energy spectrum of neutrinos in the Sun will 



study we include such kind of spin effects. We utilize Calchep [37 



De different. Thus in our 
to generate events of 



neutralino annihilation based on exact matrix elements calculation. We then import these 

□ 

events into Monte Carlo simulation program - WimpSim [39| which performs a thorough 



analysis for neutrinos propagation from the Sun to detector at the surface of the Earth 
The program contains some important features similar to those of Pythia 



40| 



4l| . Darksusy 
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42l | and nusigma 43|]. Thus WimpSim can analyze details of neutrino 



detection including both neutrino oscillations and interactions 



36 



40 
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propagation and 



45j. The detail 



Monte Carlo simulations allow us tracing every neutrino event and determining final fi signals 



with the concrete energy and angle in 
environment of IceCube detector 32, 



brmation. At the same time, we will consider the real 



46( 1 . and obtain the realistic final muon spectrum. 



This paper is organized as following. In section II we will briefiy discuss some features 
of AMSB model and the properties of the neutralino. In section III, we perform the pa- 
rameter space scan to determine different neutralino capture rates for the Sun. In section 
IV, we explore all important neutralino annihilation channels, then calculate the spectrum 
of neutrinos from the production of neutralino annihilation. In section V, we investigate 
the propagation of neutrinos from center of the Sun to the Earth. In section VI, the final 
spectrum of fi is presented. The last section contains our conclusions and discussions. 

II. NEUTRALINO IN THE AMSB MODEL 



Anomaly-mediated contributions to SUSY breaking usually appear in supergravity the- 
ory, but they are loop suppressed compared to those of gravity-mediated. However the latter 
contributions are assumed to be zero in the AMSB scenario. In this scenario there are no 
direct interaction between hidden and visible sectors [ll|, 13] • Thus the masses of neutralino 
and slepton are zero at tree level. The soft SUSY breaking terms are related to the super- 
conformal anomaly, and they appear at loop level. In the hidden sector, an auxiliary field as 



a super-gravity ground can be thought as the only origin to break SUSY 



conformal lagrangian, a compensator superfield (p is introduced 11, 
that get non zero value as 

<h=i+<F.> e\ 



12 



n order to have a 



47|. We can assume 



(2) 



then expand (j) in background value < >= (^3/2 is related to the mass of gravitino) 

and perform calculation with re-scaling the coupling of matter field. The soft SUSY breaking 
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terms related to gauginos and sleptons are calculated to be 13j: 



Ay — 

where 7, (3 are functions defined as 

d\xiZ 

7 = 



Mx = — m3/2 
9 



m 



Q 



dl . 2 



y 



(3) 
(4) 
(5) 



dg 



dy 



(6) 



dlnfj dlnjj din fx 

From these formulas we can see that the masses of SUSY particles in AMSB scenario are 
determined by ^3/2 and evolutions of gauge and Yukawa couplings. It is interesting to 



note that the mass o: 
proposed 
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epton in the Eqn.Q is negative. In literature several solutions are 



50|. Phenomenolo gica . 



terms ttt-q to the right hand of Eqn.(j4]) 13 



ly one can simply add a universal positive mass 



51| . The resulting model is dubbed as minimal 



AMSB (mAMSB) model which only depends on 



"^3/2, fno, tan/3, sign{fj,) 



(7) 



as the extra free parameters. 

The mAMSB model is highly predictive. One of the important predictions is the relation- 
ship between gaugino masses as Mi : M2 : IM3I = 2.8 : 1 : 7.1 at the low energy scale, which 
implies that the lightest neutralino is mostly wino-like over most of the MSSM parameter 
space. It is known that in MSSM the neutralino is a linear combination of gauginos and 
higgsinos, and the lightest neutralino is given as 



(8) 



On the contrary the lightest one is bino-like over most of the mSUGRA parameter space 
with GUT assumption. 

In Fig. [U we show the points which can result in the correct thermal relic density in 7713/2 
and tuq plane. We utilize Suspect code 5^ to scan parameter space in mAMSB model. For 
the case of ^"^3 + ^"^4 < 0.1, the relic density can be approximately given as |53| 

M2 ,0 



Qtkh' ~ 0.02( 



ITeV 



(9) 
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FIG. 1: The allowed points which can induce correct thermal relic density 0.085 < ^th^'^ < 0.119 
in the plane of ^-iid ^o- Curves of N^^ + N-^^ (as a measure for the lightest neutralino Higgsino 
components) are also shown. The left-top region is excluded because it can not radiatively lead to 
EWSB, and the right-bottom region is excluded because it has tachyon. The allowed points are 
generated by MicrOMEGAs [5|]- 



Such kind of lightest neutralinos usually has large annihilation cross sections in the ther- 
modynamic equilibrium of early universe, therefore it is not easy to induce the correct relic 
density. However one can still obtain the correct density if lifting the mass of neutralino 
to be as heavy as 1.9 ~ 2.3 TeV jisl . If wino-like neutralino with SU{2)l charge is 
much heavier than the weak gauge boson, the weak interaction is a long-distance force for 
non-relativistic two-bodies states of such particles. If this non-perturbative effect (namely 
Sommerfeld enhancement) of the dark matter at the freeze-out temperature is taken into 
account, the abundance can be reduced by ~ 50%. Thus the neutralino mass would be as 

n 

heavy as 2.7 ~ 3.0 TeVllTj. Such heavy neutralino may be detected through dark matter 



search experiments 15|, l55|, l56|] , but it is difficult to detect and study it at the LHC and the 
planning ILC. For the case of Nf^+Nf^ > 0.1, the light neutralino can also induce the correct 
density for a tiny region of parameter space, which lies in the vicinity of boundary of the 
electroweak symmetry breaking (EWSB). Here the decreases and Higgsino components 
become significant. 
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Now that the hght neutrahno hardly accounts for thermal relic density in the mAMSB 
model, another approach is proposed, i.e. the neutralino is not mainly from thermal pro- 
duction. Some authors discussed the mechanism that the LSP is produced by the decay of 
gravitinos Q, [l^ (or moduli fields [l^) in the AMSB model. From Eqn.Q we can see that 
neutralino is much lighter than that of gravitino. If the gravitinos have a relatively short 
lifetime, this mechanism will not destroy the success of big-bang nucleosynthesis (BBN) 1^ . 
In this scenario, the light LSP with large wino component can be the good candidate for 
dark matter. However it should be pointed out that although this kind of LSP with large 
cross section could be produced in non-thermal process, it is also constrained by cosmolog- 
ical observations. For example, in the center of our galaxy, we may observe the signals of 



;he dark matter annihilation 



57 



58 



59| due to their large cross sections. Recently, Ref. 



60( 1 reported that if assuming an excess of microwave emission in the inner Milky Way 



is due to dark matter annihilation, the synchrotron measurements give constraint to the 
annihilation rate of neutralino. The constraint depends on the dark matter distribution 
profile. More stringent constraint has been reported by Ref. Gl]- The residual annihilation 
of dark matter would produce ^Li during the BBN. Then the observations of abundance of 
^Li could constrain the dark matter annihilation rate. The low mass wino-like dark matter 
below 250 GeV in AMSB model has been ruled out 

'j le high energy collider 



wino-like neutralino at t 



experiments 



n 



64 



6l|. It is still possible to detect heavier 



62| and non-accelerator dark matter search 



65|. 



The general relation between /i and other parameters in SUSY at tree level is given by 



/i 



mjj — rri^j tan^ {3 1 



tan^ /5 - 1 



-rrir 



(10) 



As discussion in Ref. 



5l| . pure anomaly-mediated value of Tn\j^ is renormaliztion group (RG) 



invariant, and non anomaly- mediated contribution of Sm'jj^ would be even zero for some 
parameters. This focusing behavior near the weak scale makes rnjj^ only lightly dependent 
on ultra-violet (UV) boundary and would give suitable small fi for naturalness requirement. 
At the same time it will lead to considerable higgsino component of the LSP, which is similar 
as that in 'focus point region' in mSUGRA model In Fig. □ and Fig. M, 

the parameter scan was done by usin g M icrOMEGAs |54l | and Suspect 52|] with parameters 
tan/3 = 10,/i > O,mtop = ll^.QGeV (66|. In FiglH the left-top region is excluded because 
the radiative EWSB condition fails. With the decrement of the Higgsino composition 
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FIG. 2: Points which can induce rehc density lower than 0.119 in the plane of lightest neutralino 
mass and Nf^ + Nf^. Here < uiq < 20TeV and < < ^OOTeV. The curves with different 
fi at 5000, 3000, 1000, 700, 400 GeV are also plotted. We can see that the hght neutralino with 
small n is easy to get large Higgsino component. 

becomes larger near the boundary. In this parameter region, the lightest neutralino could 
be partially non-thermal-produced. At the same time, some points can even directly satisfy 
thermal relic density due to relative small wino content. In Fig. [2] we show the points which 
can induce relic density lower than 0.119 in the plane of m^o and Nf^ + Here we choose 
< rriQ < 20 TeV and < ^3/2 < 400 TeV, which is smaller than those in Fig. [H The 
reason is that the heavy neutralino, which corresponds to large mo and 7713/2, usually have 
small spin-dependent cross-section aso (see Fig. H] below )jConsequently low capture rate in 
the Sun. We impose the constraints as rrih > llAAGeV 67 1, relic density fi/i^ < 0.119 
Ao < 2 X 10-3, B{b 57) = (355±26) x 10" 



and Aaf,{Exp-SM) = (26.1±9.4) x 10' 
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69[ within 3a uncertainty. Note that in scattered plots in Fig. [3], Fig. H] and Fig. [5l we also 
choose < mo < 20TeV, < m3/2 < 4:00TeV and with the same parameters and constraints 
as in Fig. [2l In our chosen parameter space, the relic density is well under upper limit of 
WMAP. A few words on the constraint from is worth to mention. Recently, the authors 



of Ref. 



70| updated SM prediction of which has a smaller uncertainty and corresponds 



to a 3.4(7 deviation from the measured value in Ref. 69[. In mAMSB model, parameter 
space with large mo and m3/2 is heavily constrained because the SUSY contributions are 
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suppressed by heavy superparticle masses. In mSUGRA model, for the same reason the 
parameter space with large mg and mi/2 is also severely constrained. 

III. CAPTURE AND ANNIHILATION RATE OF DARK MATTER IN THE 

AMSB MODEL 

There are dark matter in our galactic halo, and the dark matter particles will scatter off 
nucleus in astrophysical objects such as the Sun and the Earth and be gravitationally trapped 
Once captured, dark matter particles will have a larger annihilation rate due 



in them 



to the larger density and will produce the energetic SM particles. These SM particles are 
mostly absorbed by the matter except neutrinos which interact feebly with ordinary matter. 
The evolution of dark matter number in the objects can be written as 

N = Cq- CaN^. (11) 

Here Cq is the capture rate and Ca =< (^v > /Veff is thermally averaged annihilation cross 
section per volume. The annihilation rate can be solved as 

r = ^CaN' = lCotanh2(v/C^t), (12) 

where t is the age of this system and | is due to the two LSPs annihilation. For the Sun, 
\/CqCa^ ^ 1 leads an approximation of above result as F = \Cq. The capture rate for the 



Sun is then given as 71 1 



\).6Lxev jcmr Viocai lU po m^o 

where piocai is the local mass density of dark matter, and viocai is velocity of these particles. 
n our numerical evaluation we adopt the values piocai = 0.3 GeV/cm^ and viocai = 220 km/ s 
Here a so and a si are spin-dependent and spin-independent scattering cross sections 
of LSP with nucleus respectively. 

In the capture process, LSPs lose their energy by scattering off the nucleus in the Sun 
and then trapped by gravity. The capture rate depends on the scattering cross section 
of LSP with nucleus. The spin-independent scattering cross section asi has been strictly 
constrained by direct detection on the Earth. For example XENON 10 set an upper limit 
for the WIMP-nucleon spin-independent cross section of 8.8 x 10~®p6 for a WIMP mass of 



lOOGel/, and 4.5 x 10" V for a WIMP mass of 30Ge1^ at 90% confidence level [Tsj. 
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The cross section for scattering of a neutralino off of a proton via spin-dependent inter- 
action is 

u,d,s 



where 



C^SD 



7i{m^i + nip) 



(14) 



m 



w 



2 

mil 

Xo 



raid? 



q q 



mi 2m 



+ [T^^Nu - tan^H/lT^L " e,)iVii]' + tan^ ^lyeJiV^, . (15) 



''W 



Here is the quark mass, dg = —Nr^/ cos/3 for down-type quarks, dg = —Nu/ sin/3 for 
up-type quarks, T^j^ is the weak isospin of the quark and Cq is quark electric charge. The 
quantity Ag measures the fraction of the nucleon spi n carrie d by the quark. The values are 



taken as Am = 0.78, M = -0.5 and As = -0.16 [7^, 175|, 176|, 177| . The first part of the Eqn. 
[T5] corresponds to the process of exchanging Z boson, and the second part of exchanging 
squark. For heavy squark, the second term can be safely neglected. 




FIG. 3: Contours of asD in the {1^^/21 ^o) plane. The constraints are the same as those in Fig. [2l 
Four benchmark points A, B, C and D, which are used in Tab. [J and Tab. |TI1 are also depicted. 

The cross sections of ctsd in AMSB model are plotted in Fig. [3] and Fig. |H In these two 
figures the largest neutralino mass is l.lTeV. For heavier neutralino the capture rate is 
rather low because asD is smaller and it is also suppressed by its mass. In Fig. HI the asn 
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FIG. 4: Allowed parameter points in the plane of lightest neutralino mass and asD- The parameters 
tan /3, sign{iJ,) , mtop are the same as those in Fig.O The constraints are the same as those in Fig.[2l 



varies from 10~^pb to 10~^pb. Note that in Fig. [3] the lightest neutralino, which lies in the 
vicinity of EWSB boundary, has higher higgsino fraction thus larger cross section aso- In 
Fig. [31 the lower left portion is excluded by direct accelerator limits on sparticle masses, 
light Higgs mass constraint. 



IV. PRODUCTION AND PROPAGATION OF NEUTRINOS 



A. Production of neutrinos in the Sun 



The dark matter annihilation in the Sun is calculated in the static limit, i.e. the dark 
matter is assumed to be at rest, because the velocity of these particles is about 220km/ s 72\ 
in the solar system. The neutrinos can be directly produced from the lightest neutralinos 
xl annihilation. However the neutrino flux is usually negligible compared to those from Xo 
annihilation to W^W~ , ZZ, Zh, f f . In this paper the direct neutrino production from Xo 
annihilation is not included. In mSUGRA model, usually contributions from each channel to 
neutrino flux play the same important role. However in AMSB model, XqS mainly annihilate 
to W^W~ which is usually 10-1000 times greater than other modes. The reason is that xh 



12 



is mostly constituted by wino and chargino is also composed nearly of charged wino, thus 
the coupling W^x'j'xf is much larger than that of ZxiXi- Note that annihilation rate of the 

1 2 1 2 1 2 1 2 

latter one is proportional to |A''i4| — |A''i3| . In Fig. [21 we have shown that |A''i4| + |A''i3| is 
small over most of the parameter space, thus the coupling ZxiXi is usually small. Especially 
in the vicinity of EWSB boundary, in which Xo can have large component of Higgsino, the 
fraction of Xo annihilation to W^W~ is always much more larger than 90%. In Fig. \^ 
the fractions of channel W~^W~ , ZZ, Zh,tt are plotted. From the figure we can see that 
the dominant annihilation channel is W^W~. For other channels ZZ, Zh,tt, only very few 
parameter points can have fraction above 10%. 
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FIG. 5: Fractions for neutralinos annihilate into WW, ZZ, ZH and tt as a function of neutralino 
mass. 

The energy distribution of induced neutrino and anti-neutrino from xh annihilation is the 
same, thus we can only consider the neutrinos. Because we plan to investigate the neutrino 
observation by IceCube which has the threshold energy around 50 GeV, we will focus on the 
flux of high energy neutrino with energy greater than 50 GeV. The production of high energy 
z/^, z/g neutrinos are the same. The reason is that the produced // is stopped by the matter 



in the Sun before its decay 781], thus it mainly produces neutrinos with low energy which 
we do not consider in this paper. Neutrinos can also be produced from charged pions from 
hadronization of partons, but they contribute mainly to low energy spectrum. Therefore for 
high energy neutrinos they are small compared with those from gauge boson and top quark 
decays. The production of Ur are usually larger than z/^, z/g because r lepton decays before 
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it loses energy 25|]. In summary the neutrino flavor ratios are given by 25 1 



Ve : Vu : Vr : Ve : Vu : Vr = 1 : 1 : r : 1 : 1 : r, 



(16) 



where r varies with the energy of neutrino. 



We utihze Calchep 



37 



to generate events, and then use Pythia 



4l| to handle the 



decays of intermediate particles contained in the events. 



1. Annihilation channel Xoxl ~^ W'^W ,ZZ 

For the annihilation channel XoXo ^ W~^W~, ZZ, the polarized effects of gauge bosons 
should be considered, which makes the energy distribution of final state neutrinos different 



from those without including polarization information [26||. In this channel, gauge bosons 
are transverse polarized with equal up- and down-type polarizations. 

In Fig. El the energy distributions of neutrino are shown for several typical neutralino 
masses. It is obvious that the polarization effects can change the shape of the energy 
distributions. In the figure the solid lines in the bottom are the analytical results of neutrinos 



which come directly from W gauge boson leptonic decay [26|. The dashed-dot lines in the 
middle are electron neutrino. We can see that the monte carlo results are higher than 
analytical results because we include the contributions from varies particle decays. Another 
distinction between analytical results and monte carlo results is that the latter ones have 
a higher tail in the low energy region because of W gauge boson hadronic decay and tau 
lepton decay. The dash lines on the top are the tau neutrinos, which are higher than those 
of electron neutrinos duo to the contributions from tau lepton decay. The muon neutrino 
has the same distribution as electron neutrino. Note that in the leptonic decay of W gauge 
boson, the muon is set to be stable in our simulations, i.e. the neutrino produced by the 
muon is not counted. The reason is that the fi loses energy quickly in the Sun and almost at 



rest before it decays into neutrino 



781 ]. The decay of tau lepton has significant contributions 



to 



;au neutrino compared with analytical result which only include neutrino from W boson 



25 



40[ |. It will enhance the number of muon neutrino via the oscillation effect [40 1. 
In Fig. U\ the energy distributions of neutrinos from annihilation channel XoXo ~^ ^^e 
plotted. For 110 GeV neutralino, the ZZ channel has narrower energy distribution compare 
to that of WW because Z is heavier than W. For heavier neutralino such different energy 
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FIG. 6: Energy distributions of neutrinos from annihilation channel XqXq ~^ W~^W with neu- 
tralino mass 110, 285 and 510 GeV respectively. For the right-top plot, the polarization effect of 



..t« t« ettetgy x = 



W gauge bosons are not included. The horizontal axis is energy x = — —. The vertical axis is 



'N~dx'^ where N is the number of annihilation XqXq ~^ W^W~ and A'^;^ is the number of neutri- 
nos. The dashed-dot lines, the short dashed lines and the dash lines represent Uf.^u^^Ur neutrinos 
respectively. The solid lines stand for the analytical results of neutrinos which come directly from 
leptonic decay of polarized W gauge boson produced in XqXq ~^ W~^W~ |26l ]. 



distribution between Z and W will be smaller, j^^^ in ZZ channel is higher than that in 
WW channel because there are two Z bosons and each Z can contribute to both neutrino 
and anti- neutrino while W boson can not. 
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FIG. 7: Same as Fig. [6] but for XoXo ~^ 
2. Annihilation channel XqXq ~^ 



In annihilation channel XqXo ~^ Z boson is longitudinal polarized 26|]. The 

energy distributions of neutrinos are presented in Fig. [HI We also compared the analytic 
results with our full simulations which include both polarization effects and various particle 
decays. Once again, the full simulations have higher tau neutrino production, right energy 
distributions shape after including polarization effects, as well as the higher low energy tail 
from lepton decay and quark hadronic decays. 

3. Annihilation channel XqXq — > tt 

For annihilation channel xJxq — ^ tt-, it has small contributions to our final results. For 
simplicity we utilize WimpSim to generate the signals with averaged t-quark polarizations 



391] . The energy distributions are shown in Fig. O 



B. The propagation of neutrino 



The neutrinos produced at the solar center will interact with the nucleus in the Sun during 
their propagation to the solar surface. The matter effects include the neutral current (NC) 
effects, the charged current (CC) effects and tau neutrino Vj. re-injection from secondary tau 
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FIG. 8: Same as Fig. [6] but for XqXo ~^ 



lepton decay. From solar center to the Earth the oscillation effects of neutrinos are also 
important and must be taken into account. In the density matrix method, the evolution 
equation is [25 1 

^ = -i\H 1 + ^ + ^ + ^ (17) 
dr ' dr dr dr 

Here the first term at the RHS describes the vacuum oscillation. The second and third 
terms describe NC and CC effect separately, including neutrino absorption, scattering and 



Ut- re-injection. The last term 



5{r)5ijj^-j^ is the initial neutrino at the production 



point. 

We use Monte Carlo program WimpSim which is event-based code to handle the propa- 
gation process of neutrinos from the Sun center to a distance of 1 astronomical unit (AU) 
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FIG. 9: Energy distributions of neutrinos from annihilation channel XqXq ~^ ^ with neutralino 
mass 285 and 510GeV respectively. The conventions are the same as those in Fig. [6] . 



391]. WimpSim treats three flavor neutrino interactions and oscillations in one framework. 



The equivalence between the Monte Carlo method and the density matrix formalism used 



m [251] has been discussed in Ref . [40|] . The detail Monte Carlo simulations allow us tracing 

every neutrino event and determining final fi signals with the concrete energy and angle 

information. Thus this approach can provide more realistic events. We take the neutrino 

□ 

oscillation parameters which given by global fit to neutrino experiments [79|] ( The 2007 
updated values in Ref. 13] have been used ). Explicitly, we choose 612 = 34.4°, ^^13 = 0° 
and 623 = 45°, the mass-squared difference Aml^ = 7.6 x lO'^e^^, Aml^ = 2.4 x lO^^eV'^ 
and CP-violating phase 6 = 0. 

V. THE DETECTION OF NEUTRINO BY ICECUBE 

When the Sun is below the horizon, the up-going muon-neutrinos may produce high 
energy muons in the Earth. When these muons travel in the ice, they will lose their energy 
via scattering with other matter. The IceCube in south pole is designed to detect such kind 
of muons. 



The analytical formula of final muon flux at the detector has been given in Ref. 26]. 
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The muon energy loss is given as 



M 

dE 



--/3£^. (18) 



where a and (3 are empirical parameters. The distance that a muon travels in the Earth 
before its energy drops below threshold energy E^^'^ , is called muon range which is given as 



R,{E) = Un{^±^). (19) 



Instead of utilizing above analytical formulas to calculate muon flux, we use WimpSim 
to get muon flux based on event-by-event Monte Carlo simulation. After obtaining the 
neutrino events at the surface of the Earth as discussed in the previous section, we then 
simulate the neutrinos propagation in the Earth and the final flux of muons. The muons 
are induced from CC interaction at actual neutrino telescope IceCube, which is mainly 
characterized by R^{E), Aejf and its geographical position. We calculate the effective 
detecting area ^e// as that in Ref. j46|. Note that ^e// in Eqn. [T]is also function of final 
muon energy. The WimpEvent uses a simple model to treat the Earth orbit and obtain the 
detector time-dependent location for each neutrino event j4o| . 

Comparing with analytical formula, the simulation can give more information about muon 
flux. The muon energy at the detector can be given by this simulation, while analytical 
formula only gives muon energy at the vertex where neutrinos convert into muons. The 
information about angle between muon velocity and the direction from the Sun to the 
detector can also be obtained which is important to reduce background. 

We show our final results for some benchmark points in Fig. [10] after multiplying initial 
dark matter annihilation rate in the core of the Sun. Here we have summed the fi^ events 
together and taken E^'^^ = 50GeV into account. After produced from charged current 
scattering process, the muon may have multiple Coulomb scattering in the ice 36|, ISjJ . Thus 
the final direction of muon velocity is not the exactly same direction from the Sun to the 
detector. From Fig. [101 "we can see the deviation of muon direction is not large, and most of 
muon events are in the scope of ±2° around the direction from the Sun to the detector. And 
the signals from high energy neutrinos produced by heavy massive neutralinos annihilation 
are less affected by Coulomb scattering as naive expectation. Note that the IceCube can 
reach angle resolution as small as 1° [32|, our results suggest a possible method to reduce 
the background events from other sources. 
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FIG. 10: The muon events (/i"^ + /U~) distributions in IceCube as a function of angle 9 which 
is defined as angle between the muon velocity and the direction from Sun to the detector. The 
IceCube 1*^ angle resolution is not reflected here, and every bin is taken as 0.2*^. The threshold 
energy for muon is taken as 50GeV. 

Now we switch to the discussions on the backgrounds. The largest backgrounds come 
from comic ray. When cosmic ray interacts in the atmosphere around the Earth, high energy 



neutrinos would be created and 
can produce muons 



36, 



82 



;ravel to the detector. During the propagation, the neutrinos 
|. Such kind of high energy neutrinos do not have special 
direction. Therefore if we only observe the events around the special direction from Sun to 
the detector, the backgrounds would be dramatically reduced. Other possible background 
are muons from cosmic ray interacting with the atmosphere around the Earth 8J] and 
neutrinos from cosmic ray interacting with particles in the Sun's corona 851]. However these 
two kinds of backgrounds are not important in the detection of muons at IceCube. For 
the rough estimate, we only take the atmosphere neutrino backgrounds into account. The 
atmosphere neutrino flux are taken from Ref. 82] and method to estimate the muon rate 
from Ref. 26|. The resulting background muons are not the function of i^-^'" but for E^, 



i.e. energy of the induced muon. It should be noted that more accurate results need Monte 
Carlo simulation. 

Our final results are depicted in Fig. [11], Tab. [T] and [Tll The four parameter points are 
almost dominant by annihilation channel WW which is one of the characters of AMSB model. 
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FIG. 11: Differential muon rate in IceCube as a function of energy of muon for several neutralino 
mass as 110, 285, 510 and 1000 GeV. We applied different angular cuts on muon zenith angle in 
order to get better significance according to the result of Fig. [lOl The red dashed (the blue dotted 
) lines represent the muon signals (the atmosphere backgrounds). 

All four benchmark points, which satisfy all the constraints from other experiments, have 
the same order of magnitude of (Jsdi which is essential for the muon detection at IceCube. 
Point A with llOGeV neutralino has much smaller signal than that of Point B with 2'$)bGeV 
neutralino because the former is affected by the energy threshold and muon range i?^. Point 
C and D with heavier neutralino are less affected by the energy threshold but the capture 
rate Cq is suppressed. We applied different angular cut on zenith angle of muon to get better 
statistical significance according to the results in Fig. [101 From Tab. [Tll we can see the signal 
annual events can reach 102 for Point B and the statistical significance astat can reach more 
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Point 


m^{GeV) 


7Tio(TeV) 


7713/2 (TeV) 


/i(GeV) 


Rh 


crsoipb) 




A 


110 


2.9 


34 


420 


4.4 % 


5.1x10"^ 


8.3x10^9 


B 


285 


6.1 


88 


498 


7.1 % 


4.2x10"^ 


1.0x10^9 


C 


510 


9.5 


162 


620 


19.2 % 


5.0x10"^ 


4.1x10^^ 


D 


1000 


15.3 


314 


1120 


14.5 % 


1.1x10"^ 


2.2x1027 



TABLE I: Several benchmark points (point A-D) in mAMSB model. Rh is the higgsino fraction 

1 2 1 2 

I -^13 1 + 1-^141 5 o'sD is the spin-dependent cross section and Cq is the neutralino capture rate for 
the Sun. 



Point 


a-ww 

O'tot 


O'ZZ 
O'tot 


O'tot 


O'tot 


ecut (°) 


Sig 


BG 


O'stat 


A 


100.0% 


0.0% 


0.0% 


0.0% 


2.0 


23.3 


24.1 


4.7 


B 


99.1% 


0.0% 


0.9% 


0.0% 


2.0 


102.1 


24.1 


20.8 


C 


95.6% 


0.2% 


4.2% 


0.0% 


1.5 


70.3 


13.6 


19.1 


D 


96.9% 


0.1% 


3.0% 


0.0% 


1.0 


3.0 


6.0 


1.2 



TABLE IL One year total muon rate at IceCube for various benchmark points (point A-D) where 
muons satisfy Ethr < < 300Gey with Ethr = 50GeV . Here ^ are the fractions 

of annihilation channels WW, ZZ, Zh, tt separately. 6cuT is the angular cut for zenith angle of 
muon. Sig is the annual muon rate and BG is the annual muon background from atmosphere 
neutrino. We define statistical significance as agtat = S/\fB. 

than 20. Based on our numerical estimation, we can obtain the rule of thumb for IceCube 
to discover high energy neutrino from Sun in AMSB model, i.e. aso > lO'^pb. It should 
be noted that Nf^ + Nf^ > 4% approximately corresponds to the condition aso > I0~^pb 
based on the results in Fig. [1] and [3l 

VI. CONCLUSIONS AND DISCUSSIONS 

In this paper in the mAMSB model we have investigated the muon event rate at IceCube. 
Such kind of muon events are induced by high energy neutrino flux from neutralino anni- 
hilation in the Sun. We studied the detail energy and angular spectrum of the final muons 
at the detector based on event-by-event simulation by using Monte Carlo code WimpSim. 
More precisely we simulated the processes since the production of neutrino via neutralino 
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annihilation in the core of the Sun, neutrino propagation from the Sun to the Earth, as well 
as the converting processes from neutrino to muon. Our results show that in the mAMSB 
model it is possible to observe the energetic muons at IceCube, provided that the lightest 
neutralio has relatively large higgsino component, as a rule of thumb Nf^ + > 4% or 
equivalently aso > 10~^pb. Especially for our benchmark Point B (see Tab. [I]), the signal 
annual events can reach 102 and the statistical significance can reach more than 20. This 
parameter space is very similar to the 'focus point' region in mSUGRA model. 

It should be emphasized that our study has much more improvement compared to pre- 
vious investigations. For example, we include both polarized effect of gauge bosons and 
secondary neutrinos from the decay of leptons and quarks. The polarization effects of the 
gauge bosons can have important influences on the neutrino energy spectrum, as shown 
in our results and in Ref. 



for example in Ref. 



261 ]. The secondary neutrinos are usually handled by Pythia, 



25| and WimpSim 39|]. But none of them include both effects. We 



use Monte Carlo code WimpSim to handle the interactions includes NC effect, CC effect 
and tau neutrino z/^- re-injection from secondary tau lepton decay. Comparing with the 
previous analytical result, the Monte Carlo approach can not only give the same evolution 
results as those by solving Eqn. [T7] but also the angular distribution of muons as shown 
in Fig. [TUJ With the angular distributions at hand, we can choose 2° as normal cut in 
range lOOGeV < m^i < 500GeV to reduce the background from atmosphere neutrinos. 
With heavier neutralinos, smaller cut can be chosen to further reduce the background. As 



shown in Ref. 



60|], the light neutralino in this scenario is disfavored. Thus the techniques 



to suppress background in order to extract low signal events for heavy neutralino is very 
important. 

Last but not least, we want to emphasize that the final energy spectra of muons will 
have similar shape as those in Fig. [11] in the mAMSB model, because the annihilation of 
neutralinos are almost XoXo ~^ W~^W~ . Thus it is possible to distinguish among AMSB 
model from other SUSY breaking scenarios due to different annihilation modes if enough 
muons can be collected. 
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